Toy Model for Pion Production in Nucleon-Nucleon Collisions 
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' I. INTRODUCTION 

G\ ' 

Interest in studies of pion production in nucleon-nucleon collisions at energies near threshold has been re-vitalized 
by the appearance of excellent high quality data Q| . The fact that low- and medium-energy strong interactions are 
controlled by chiral symmetry led to an early hope that chiral effective theories could be used to analyze these processes 
and achieve a fundamental understanding of the production process. Indeed, there are now tree level calculations 
and even loop calculations PHlC| available in the literature p| . The early excitement was quickly abated by 
the realization that proper evaluation involves surmounting several severe difficulties, which are caused by the high 
momentum transfer nature of this threshold process. The initial relative momentum between the two nucleons must 
be at least pi = ^Jm^M^ . This means that the chiral expansion is in terms of powers of y/ 111^ /Mn instead of m^/MN 
[§,[l2| , which complicates carrying out the expansion and verifying its convergence. However, issues of convergence are 
not the focus of the present work. Instead, we address some technical questions which arise during the evaluation of 
the relevant matrix elements. 

It is worthwhile to discuss some general features of the pion production process before describing our specific 
technical issues. Pion production occurs when the mutual interactions between two nucleons cause a real pion to be 
emitted. The leading term is one in which the initial and final state two-nucleon (NN) scattering allow a pion to be 
emitted by a single nucleon emission. The next tree level contribution occurs when a virtual pion of four-momentum 
q produced by one nucleon is knocked on to its mass shell by an interaction with the second nucleon. This is the 
so-called re-scattering diagram. This process typically occurs accompanied by low- momentum-transfer initial and/or 
final state interactions. The evaluation of these diagrams, including the case when the pion exchanged between the 
two nucleons may be on shell, is our focus. Our strategy will be to introduce a toy model, which is simple enough 
to allow the exact evaluation of certain amplitudes. Then we may assess various approximations by comparing the 
resulting amplitudes to the exact results. 

In general, one could obtain the necessary transition matrix elements by evaluating the relevant Feynman diagrams. 
However, the initial and final state interactions are accurately treated using an appropriate NN potential within a 
three-dimensional Schrodinger equation formulation. Thus one needs to obtain a three-dimensional formulation from 
the more general Feynman procedure. This has been done in an ad hoc manner in Refs. |^,^|||-|^]: one guesses the 
energy dependence of the virtual pion-nucleon (irN) interaction, and uses a Klein-Gordon propagator for the pion 
propagator. However, there is a general method to derive a three-dimensional theory which is equivalent to the 
Feynman diagram approach, namely the method of considering all the time-ordered diagrams — the use of time- 
ordered perturbation theory (TOPT). In this formulation, one finds only NN, irNN and ttttNN propagators in the 
tree-level re-scattering diagrams. The Feynman Klein-Gordon pion propagator does not appear explicitly. Thus our 
first focus is the appropriate propagator. In particular, we will compare different prescriptions used in the literature 
with the exact result derived in the toy model. 

Another issue to be addressed is that of the proper choice of the energy variable go of the exchanged pion. The 
value of qo is critical because the chiral 7riV interaction includes seagull vertices involving Oqtt, such as the isovector 
Weinberg- Tomozawa interaction N^tN ■ (ir x doir) and the isoscalar N^ N (doir) 2 . In case of the isoscalar re-scattering, 
which is most relevant for threshold ir° production, this seagull term is oc q^m^. Its actual size is crucial: for on-shell 
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irN scattering at threshold (go = m w ) there is an almost complete cancelation of different, individually large terms 
leading to a very small irN isoscalar scattering length E| . If one moves away from the threshold or the on-shell irN 
kinematics, however, this cancelation gets less and less effective. Thus the numerical value of the isoscalar re-scattering 
term is very sensitive to the details of the individual terms. Note that, because of the Weinberg- Tomozawa term, the 
proper choice for q is also relevant for the isovector re-scattering that contributes to charged pion production. 

If one simply evaluates the re-scattering diagram at threshold, neglecting initial- and final-state interactions, it is 
clear that go = Keeping this value fixed also when including the distortions leads to an amplitude that is 

opposite in sign to the on-shell scattering amplitude, and interferes destructively with the single-nucleon emission 
term [^]|| . This choice for go in combination with the use of the Klein-Gordon propagator for the pion will be called 
fixed kinematics approximation in what follows. Refs. found that the computed cross sections fell well below the 
data, unless many other even less- well constrained terms are included q|. However, once the nucleons are no longer 
on shell, there are other prescriptions in the literature for choosing q . In a Feynman diagram this is the difference 
between the Oth component of the nucleon four-momentum before and after pion emission. Thus one might find it 
natural to set qo equal to this difference in energies. Using this energy difference prescription in the distorted wave 
Born approximation (DWBA) calculation of pion production leads to a re-scattering diagram which also has a sign 
opposite to that of the single nucleon term, but which is about three times larger in magnitude As a result, one 
can reproduce the magnitude of the total cross section using only the re-scattering diagram. This energy prescription 
will be called (E — E') approximation below. 

In addition, having the toy model at hand, we also want to study the importance of terms that go beyond the 
DWBA, namely the so-called stretched boxes (c.f. Fig|l], diagrams F3 and F4). These necessarily occur in the 
three-dimensional framework and represent diagrams where there is no two nucleon cut. 

Note that the questions under investigation affect not only chiral perturbation theory calculations, but also more 
phenomenological approaches. For example, Ref. |l4| used the (E — E') prescription for the pion re-scattering when 
investigating the influence of nucleon resonances on the production process. In the so-called Julich model |l5[ the 
full TOPT propagator was used, but with its energy fixed to the production threshold. Thus, a clarification of these 
formal issues is necessary before one can draw conclusions about the physics of the process. This paper is meant to 
be a step in that direction. 

It is important to realize that one can not resolve the ambiguity in the choice of go or the proper nNN propagator 
(in what follows this quantity will sometimes in a somewhat sloppy way be called "pion propagator" ) by appealing to 
data. These are questions about the theory which arise due to the manner the DWBA procedure was implemented 
. Furthermore, the slow convergence of the momentum expansion requires one to resolve these difficulties before 
evaluating loop diagrams. 

One needs to construct an ab initio theory of pion production. Doing this for the realistic case requires that one 
considers several important features including i) the spin and isospin of the two-nucleon system; ii) the Goldstonc 
boson nature of the pion as an odd parity system degenerate with the vacuum; and iii) a realistic NN potential. 
However, none of these features affects directly the questions which we want to examine. Therefore it is appropriate 
to construct a toy model which is simple enough to evaluate so that exact answers can be obtained. Then we can 
consider the various choices for go & n d for the "pion propagators" as testable approximations. In the section II we 
formulate our toy model, and examine the various approximations for final- and initial-state interactions in sections 
III and IV, respectively. Our conclusions are summarized in section V. 



II. THE TOY MODEL 



The first step is to construct the necessary solvable model. Therefore: 

i) We consider the production of a scalar "pion" field which has a Yukawa coupling with the nucleons. (We shall 
leave out the quotes around pion in the following text.) 

ii) We include two nucleon fields, or, alternatively, treat nucleons as distinguishable. As a consequence, we need 
only include pion emission from one nucleon, but not the symmetric term where the pion is emitted from the other 
nucleon. We do not have to worry about several spin-isospin channels, and respective projections. The simplicity of 
the model is retained by allowing the pion to couple to only one nucleon field. As a result, the effects of pion exchange 
between two nucleons does not enter. 

iii) A focus of the paper is the pion re-scattering by one nucleon. This pion re-scattering is described by a irN 
seagull vertex which is inspired by the chiral ttN interaction Lagrangian. 

iv) In order to mock up the nuclear interactions we include the exchange of a scalar sigma field, which also couples 
to nucleons via Yukawa coupling. Since the magnitude of this coupling has nothing to do with the way to treat the 
pion energy, we consider the case of small coupling, and therefore need to only consider one sigma exchange. 
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v) Because pi/M^ = ^Jm^jM^ < 1, it is typical to treat this problem using a non-relativistic expansion. In the 
following we will examine only the leading terms in this expansion. In particular, contributions from anti-nucleons 
are not considered. 

Therefore, we consider the following toy model defined by the Lagrangian: 

C = NHid + ^)N + I [OV) 2 - ml^ + (d,a) 2 - m 2 a 2 ] 

i=l,2 N 

J* 1=1,2 ^ i=l,2 

Here Mjv is chosen as the physical nucleon mass of 939 MeV and similarly m n is taken as 139 MeV. The mass of the 
a meson and the cutoff A on the momentum integrals are taken as parameters in the theory, to be specified below. 

It is important to immediately display some of the non-realistic features of this toy model. For simplicity, we did 
not enforce chiral symmetry, which would have required a derivative coupling of the pion to nucleon spin, instead of 
the simpler Yukawa coupling. We are concerned with near-threshold kinematics so that a scalar particle is produced 
in an S wave, as is the final NN pair. Angular momentum conservation requires that the initial NN pair also be 
in an S wave. In the real world, however, pions are pseudoscalar and thus the production of S'-wave pions calls for 
a P wave in the initial state. Furthermore, the toy model includes no strong short-range repulsive NN interactions 
which keep the nucleons apart. Thus the nucleons have stronger overlap for our toy model than in a more realistic 
treatment. However, to a given order in the coupling constants we can obtain exact amplitudes for this model, and are 
therefore able to study the various treatments of qq and the irNN propagator to determine which, if any, reproduce 
the exact model answers. 

In a DWBA calculation of threshold pion production, the tree-level re-scattering diagram is influenced substantially 
by the contributions from the initial and final state interactions. In this toy model calculation we will therefore for 
simplicity concentrate on the DWBA terms where we have only initial or final state NN interactions. We will in 
this paper ignore the re-scattering diagram with DWBA contributions in both initial and final NN interactions since 
this is a two-loop integration term. Again for simplicity we will, as discussed, simulate the NN interactions with a 
single a exchange between the nucleons which occurs before or after the pion re-scattering process — the initial-state 
interaction and the final-state interaction, respectively. We will discuss these two cases separately below. In addition 
there are graphs in which a a is exchanged in between the emission and re-scattering of the virtual pion. We ignore 
these here, as they are not relevant for the issue at hand. All of our diagrams are evaluated at order j^-g 2 -^. In the 
following we do not display these factors, as well as other constants which are common to all of the amplitudes. 



III. FINAL-STATE INTERACTION 

The exchange of a a meson in the final state is given by the Feynman graph FO in Fig. [lj. We consider threshold 
kinematics in the center-of-mass frame, and use the following notation. E(E') represents the energy of a nucleon in 
the initial (final) state, with Etot = IE = 2E' + m v = m T (at threshold: E' = 0). In addition, u> q = yj m 2 , + q 2 

and uj a = \l to 2 + k 2 denote the n and a meson on-shell energies, and E" = k 2 /2Mn the energy of an intermediate 

nucleon. Here k = p + q, where p is the initial nucleon three- momentum. We choose the pion momentum q to be the 
integration variable, so that the diagram shown in Fig. |l| (FO) corresponds to the following four-dimensional integral: 

f JN_ { i 

J (2tt) 4<Z ° \(E + q -m w -E» + ie)(E" + q - E - ie) 

1 

(qo - Lu q + ie)(qo + uj q - ie) 
1 

(qo - E + E> + Lu a - ie){q - E + E' - uj a + ie) 
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FIG. 1. The diagrams that occur when the sigma exchange appears as final-state interaction. The analog diagrams 11-14 with 
the sigma exchange in the initial state are considered as well. For the first two diagrams the two possible time orderings for the 
sigma exchange lead to identical expressions. 



All DWBA calculations are made using a formalism in which matrix elements are given as three-dimensional 
integrals. Thus the first step is to find the appropriate three-dimensional expression by performing the qq integration. 
Obviously, Eq. (|J) contains three poles in the upper half plane as well as three in the lower half plane. One way to 
proceed would be to close the contour on one of the half planes and pick each of the three poles enclosed. However, it 
is more convenient to perform a partial decomposition, in which the poles of the pion propagator are isolated before 
the go integration is carried out. It should be emphasized, however, that the final result does not depend on the 
method of its evaluation. It should not come as a surprise that the final result of the qo integration agrees exactly to 
the one of TOPT, as the equivalence between the Feynman prescription and TOPT is well known. This is illustrated 
in Fig. [l], and the resulting amplitude is given by 



d 3 



q u„ 



(2tt) 3 AugLJa I {Etot - E> - m n - E» - uj a )(E t ot - m n - 2E")(E tot - E - E" - uj q ) 

2 

~ {Etot -E'-m n - E" - Ua)(Etot - m„ - 2E")(E tot - E - E v - - uj q ) 

1 

+ (E tot ~E>- m ff - E" - io a ){E tot - E> - E - uj a - uj q )(E tot - E - E" - oJ q ) 

1 

(Etot — E' -m w - E" - uj a )(E to t - E - E' - m v - cu q - Lu a )(E to t — E — E" — m w - uj q ) 



(3) 



in which the successive four terms can be immediately matched to the diagrams F1-F4. In particular, the last 
two terms are those of the stretched box diagrams which have not yet been considered in any calculation for pion 
production. We will examine their importance below. Note, that there is no freedom with respect to what is the 
appropriate choice for qo in the numerator of Eq. (||) . The pole structure of Eq. (||) in combination with the way the 
partial decomposition was performed forces qo = uj q in Eq. (|3|), which is an exact equation. 

To compare Eq. (|^) to expressions used in the literature it is useful to combine the first two lines to obtain the 
final state interaction contribution to the DWBA amplitude: 



I 



d 3 q 

(2tt) 3 '° V Etot -m^~ 2E" J 2 



k_qTOPT 



7T 5 



(4) 



where the sigma potential is 



V a (k 2 ) 



u a (E t ot — E' — m n — E" - uj a ) ' 
and the TOPT irNN propagator — the exact propagator — is given by 

f 



G 



TOPT 



>2 / fc2 

1 ^1 + 2U^ 



(5) 



(6) 



Apart from the k 2 term in the TOPT irNN propagator, Eq. (|]) agrees with what is known as fixed kinematics 
approximation PJq]. As was explained above, this approximation is defined by the use of 772^/2 for the pion energy 
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in both in the ttN seagull vertex and in the pionic Klein-Gordon propagator. In the realistic case (when appropriate 
nucleon wave functions are used for the distortions) the significance of the k 2 /2Mjy term in the pion propagator of Eq. 
(||) can be estimated by noting that k 2 /2M^ is of the order of the the off-shellness of the intermediate nucleons. Since 
the final state is at rest we can estimate k 2 = 0(m 2 ) [{L2|| . It then follows in the absence of initial-state interactions 
that the loop three-momentum is \q\ ~ Pi. We can expand Eq. (||) in powers of m^/M^, and get 



Here the Klein-Gordon propagator in the fixed kinematics approximation is defined by 

G« G = \ • (8) 

The right-hand side of Eq. is already expressed in terms of the expansion parameter of the underlying effective field 

theory, Thus —at the level of accuracy accessible today — we expect this Klein-Gordon propagator to 

be a good approximation for those diagrams where the NN interaction appears in the final state. Such considerations 
are not necessarily germane here, however, as we have not enforced the chiral symmetry on which power counting is 
based. The physical scales appearing in the final state of this model are set by the parameters A and m a , which we 
take to vary over a large range. 

Let us now discuss the numerical significance of the individual terms above. Our toy model allows us to answer the 
following three questions: 

• What is the relative importance of the stretched boxes (F3 and F4 in Fig. |l|) compared to the "DWBA- 
contributions" (Fl and F2)? 

• How good an approximation is the propagator G^ G of Eq. (||) compared to the exact propagator G™ PT of 
Eq. (§)? 

• What is the effect of different treatments of the pion energy "go" at the irN seagull vertex (fixed kinematics 
compared to the (E — E') prescription)? 
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FIG. 2. Importance of the stretched boxes for different choices of the cutoff as a function of the mass of the a meson. The 
ratios of the stretched boxes with respect to the DWBA piece, Eq. aA), are shown for A — 3ra„- (solid line), Wm^ (dashed line) 
and oo (dot-dashed line). 
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The answer to the first question is obviously a function of the a mass, since the DWBA contributions should lead 
to results that are proportional to (pi/m a ) 2 , whereas the stretched boxes lead to (pi/m^) 4 . In Fig. || we show the 
ratio of the stretched box contributions to the DWBA part as a function of the mass of the sigma meson. The three 
curves correspond to three different values of the cutoff for the radial integration. As expected, the curves all fall as 
for large m a . The strength of the stretched boxes never exceeds 6%. This justifies a DWBA treatment of the 
final state in this pion production process. 

The answer to the second question is presented in the left panel of Fig. ^ as a function of the cutoff in the momentum 
integration. We evaluated the DWBA piece, Eq. (Q), with the exact propagator (^) and with the approximate 
propagator (||). The mass of the sigma was chosen to be m a = 550 MeV. The solid curve shows the result using the 
approximate pion propagator in units of the exact result. Thus, the deviation of this approximation from the exact 
result never exceeds 25%. 

FSI only ISI only 




02468 10 02468 10 

A/m n (A-p,)/m, 

FIG. 3. Effects of the different approximations to the "production operator" for different cutoffs. The left (right) panel shows 
the result for a sigma exchange in the final (initial) state in units of the exact answer of Eq. (U) (Eq. \l ]y). Using the "fixed 
kinematics approximation" pion propagator leads to the solid curve. The dashed curve is the result when using the (E — E') 
approximation. 

The third question goes to the choice of energy variable "go" at the nN seagull vertex reported in the literature [Q. 
To simulate this choice we replace "go" = m^/2 in the numerator of Eq.(|]) with 

«q »=E-E», (9) 

together with the approximate pion Klein-Gordon propagator. This was defined as the (E — E') prescription above. 
In many reactions Eq. (^) is an appropriate replacement, because the nucleons remain almost on the mass shell in 
the intermediate states. However, as soon as large intermediate momenta are accessible, this treatment might be 
questionable. In Fig. || the dashed curve shows the results of using the (E — E') prescription, again in units of the 
exact result. Within our toy model, the result shows that this approximation is not reasonable for calculations of 
threshold pion production. Note that this amplitude is very sensitive to the sigma mass, which acts as a regulator. 
If the sigma mass is taken to be larger, then the result changes even more dramatically with the cutoff. A change in 
sign happens at the point where the cutoff is big enough for the effect of the k 2 /2M to overcome that of m^/2 (the 
larger the intermediate momentum, the larger E" ) . 

The net result of the toy model for the final-state interaction case is that using "go" = m^/2 in both the virtual 
ttN seagull scattering vertex (numerator) and in the approximate pion Klein-Gordon propagator is very reasonable. 

IV. INITIAL-STATE INTERACTION 

We now consider the case when the sigma exchange occurs before the re-scattering process. A reduction to the 
three-dimensional integral (or starting with the TOPT expression) gives the following four terms: 



(i 



M = [ d3q UJq 

1 J (2tt) 3 4u q tJ a \ (Etot -E' — E" - uj q ){E tot - 2E")(E tot -E-E*-w 

2 



(E tot -E>- E" -u q - m^){E tot - 2E")(E tot — E — E" — w a ) 

1 

' {Etot -E>- E" - w q )(E tot -E-E'-u q - w a )(E tot - E - E" - cu a ) 

1 

" (Etot — E' — E" -Lj q - m*){E tat — E — E' — uj q — uj a — m^(E tot — E — E" — uj a ) 



, (10) 



where again oj q = \/m\ + q 2 ,u) a — yfftf + (<f +p) 2 , and E' = 0, but the energy of an intermediate nucleon is 



N- 



E" = q 2 /2M,, 

As before, the first two terms in Eq. ( |l0| ) correspond to box diagrams and the last two to stretched boxes, cf. Fig. 
|l|. In case of the initial-state interaction the stretched boxes still turn out to be smaller than the boxes, but less so: 
~ 30%. Note that this is also of the size expected in the real world where the expansion parameter of the EFT is 
v/]§^ — 0.4 PJl^|. We therefore concentrate on those terms containing the NN propagator, Gnn — [Etot — 2E") 
only and obtain 



-l 



d 3 q 



Gl% T ( = ' ^ ) V„ (11) 



where 



(27r)3 2 ™ v \Etot-2E" 



K(g 2 ,fc 2 )= -, r, (12) 



and in the initial-state interaction case the TOPT irNN propagator reads 

topt 1 

- ' (13) 

This looks like a DWBA expression using the prescription. Due to the large initial momentum, the unitarity 

cut of Gnn turns out to be an essential feature. 

Similar to the section on final-state interaction, we investigate the fixed kinematics approximation and the E — E' 
approximation using the free pionic Klein-Gordon propagator, Eq. (g). This means especially that in the G™^ T of 
Eq.(|l3"|) we set E" — 771^/2, which implies on-shell intermediate nucleons: q 2 = m^MN- In the E — E' approximation 
we further replace "go" = m v /2 by E — £"'. In the right panel of Fig. || we demonstrate the inadequacy of both 
approximations, compared to the exact result given by Eq. (|TT|). 

Due to the large initial momentum the imaginary part of these diagrams turns out to be of the order of the real 
part. (Since we work at the kinematical threshold of pion production the imaginary part from G^nn is zero). Since 
all the approximations where constructed such that they agree once the intermediate two-nucleon state goes on-shell, 
all the individual results agree for the imaginary part. 

The question becomes why do both approximations show such a large deviation from the exact result of Eq. (|l] 
The cause can be traced back to the appearance of a TtNN cut in the exact propagator: from Eq. (|l3|) we see that the 
propagator G™^ T diverges as \q\~ 2 when q approaches 0. On the other hand for small q the free pionic Klein-Gordon 
propagator G% G goes to a constant. It is the very different nature of the infrared behaviors of the propagators, G^ G 
and G™^ T , that leads to the large deviation of (the real part of) the amplitude from the result of Eq. (pi]). 

Having identified the nNN cut as an important feature of the production reaction, a natural question that arises is 
how to set up a counting scheme capable of covering this. Note that contrary to the more conventional contributions, 
where the scale of typical momenta is set by the initial momentum pi = y/MNTn^, the irNN cut pronounces momenta 
of the order of the external pion momentum. It can not be a part of a toy-model investigation to completely resolve 
this matter — after all our model interaction is not consistent with the requirements of chiral symmetry. However, we 
will use the last part of this section to suggest a possible method to address the issue. 

To this end we will rewrite Eq. ([n]) such that we isolate both the NN and the wNN singularities. For this purpose 
we are guided by the unitarity transformation method of Ref . ]l6| , |l7| ] . This method is one way to isolate the different 
singularities of a particular diagram. In this case the scattering amplitude can be written as (where for clarity we 
suppress J d 3 q as well as some over all factors) [pf: 
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Mi = M NN + M %NN + ■■■, (14) 
■"v.v l — ( -^-^ 1 ■ (ir,-. 




where the ellipsis denote the stretched box TOPT diagram contributions. We see that the above amplitude has two 
physical singularities due to the NN and irNN scattering states. We find numerically that M^nn is about 5 times 
larger than Mnn when evaluated with a cutoff of A = 10m w . This large effect of the ttNN cut in the toy model is 
also responsible for the stronger effect of initial-state stretched boxes, as the latter also contains the irNN cut (see 
fourth line of Eq. ([h])). These two points highlight the numerical significance of the three-particle cut. 

Note that this is not a unique separation of the two branch cuts. To make closer contact with previous work [ p~6 |l7 
we can rewrite Eq. ( fTi]) in a form closer to the (E — E') prescription: 



Mj = M' NN + M' 7rNN + • • • , (17) 

1 / E" 
-2E" \ E ,,2 ~uj 2 



M ™ = V., (18) 



m n — E — u) q \ h i A — ui~ 



Clearly Mnn + M^nn = M' NN + A1^ NN , although some shift of strength is then achieved between TV TV and irNN 
contributions: in this case we find the contribution from M^ NN larger in magnitude than M' NN by a factor 2, using 
the same cutoff A = lOm^. It remains to be seen which splitting is the most appropriate in the realistic case. 

The most significant finding for the case of the initial-state interaction is therefore that in the toy model the three- 
body ttNN branch cut of G^nn is very important. The importance of this cut has been advocated before, for example 
in Ref. [ fl9| |. Here the static propagator, which was defined as being part of the fixed kinematics approximation as well 
as of the (E — E') approximation, leads to erroneous results for the real part of the amplitude. 

However, it is important to remark that we expect the importance of this branch cut to be much smaller in the 
real world. Indeed, as we have seen, close to threshold this type of contribution comes from three-momenta near 
0. In the real world, chiral symmetry suppresses such contributions. The pion coupling in leading order in chiral 
perturbation theory, for example, goes through the pion three-momentum. In our toy, chiral symmetry does not play 
a role, the pion coupling is a simple Yukawa coupling, and both initial and final NN states are in relative S waves, 
which enhances the influence of the irNN cut. The power counting developed in Refs. |2 12 does take into account 



chiral symmetry — thus the correct factors of momenta — and suggests a suppression of these branch effects, as long 
as the momentum of the emitted pion is 0{m 7T ) (or less). Clearly, it is important to further study the power counting 
— in particular in conjunction with the unitary transformation method — in the realistic case. 



V. SUMMARY AND CONCLUSION 



We have investigated various approximations for pion production by defining a toy model which allows the com- 
putation of exact model transition-matrix elements. Because it lacks chiral symmetry, this model has the unrealistic 
features that both the initial and final state JVJV wave functions are S states. The influence of NN correlations 
that suppress the short-distance wave functions are absent from the toy model. Furthermore, diagrams with both 
initial- and final-state interactions could also be important in more realistic calculations. We have performed some 
test calculations using the Reid NN potential, which indicate that NN correlations do modify some of the toy model 
findings at a quantitative level. However, the toy model allows the compilation of exact results at a given order in 
the couplings, and thus some qualitative insight. 

The findings of this paper can be summarized as follows: 

• The stretched box contributions are numerically small compared with boxes. 

• For the final-state interaction, only the fixed kinematics approximation (for both propagator and vertex) turns 
out to be appropriate. 
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• If a loop with the initial-state interaction is included, the contribution of the nNN cut is very important and 
has to be taken into account properly, which is not done in the common approximations. 

The first two findings are in accord with the expectation from the existing power counting for pion production in 
the effective field theory PjKq]- Indeed, according to this power counting stretched boxes involving pions are sub- 
leading and those involving heavier mesons are absorbed in higher-order local operators. Moreover, due to infrared 
enhancements that lead to the (quasi) bound state in the NN interaction, the effect of the final-state interaction in 
realistic calculations should be by far dominant close to threshold. 

The third finding is perhaps surprising. However, chiral symmetry is expected to be crucial in suppressing this 
contribution in the real world, because the ttNN cut emphasizes small momenta. Clearly, the importance of the 
three-body nature of the intermediate state needs to be further examined in realistic calculations. 
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